RECTANGULAR i?-TRANSFORM AS THE LIMIT OF 
RECTANGULAR SPHERICAL INTEGRALS 



FLORENT BENAYCH-GEORGES 



Abstract. In this paper, we connect rectangular free probability theory and spherical 
integrals. We prove the analogue, for rectangular or square non-Hermitian matrices, 
of a result that Guionnet and Maida proved for Hermitian matrices in [GMOSj . More 
specifically, we study the limit, as n, m tend to infinity, of — log E{exp [y' nmOXn]} , where 
G M, Xn is the real part of an entry of UnMnVm, Mn is a certain n x m deterministic 
matrix and Vm are independent Haar-distributed orthogonal or unitary matrices with 
respective sizes n x m x m. We prove that when the singular law of M„ converges 
to a probability measure /x, for 6 small enough, this limit actually exists and can be 
expressed with the rectangular i?-transform of This gives an interpretation of this 
transform, which linearizes the rectangular free convolution, as the limit of a sequence 
of log-Laplace transforms. 

Introduction 

In this article, we study the hmit, as n, m tend to infinity in such a way that n/m tends 
to a hmit A G [0, 1], of 

- logE{exp[v/^e 3f?(Tr(E„f/„M„Kn))]}, 
n 

where ^ e R, Mn is a certain n x m deterministic matrix, Un, Vm are independent Haar- 
distributed orthogonal or unitary matrices with respective sizes n x n, m x m, En is an 
m X n elementary matrix (i.e. a matrix which entries are all zero, except one of them, 
which is equal to one) and 3fJ(-) denotes the real part. 

The departure point of this study is the work of Collins, Zinn- Justin, Zuber, Guionnet, 
Maida, Sniady, Mingo and Speicher who proved, in the papers |ZZ03l IC031 IGM051 ICS061 
ICS071 ICMSS07] , that under various hypotheses on some n x n matrices An and Bn and 
a positive exponent a, the asymptotics of 

^ log E{exp[ne Tr(5„f/„A„f/:)]} 

are related to free probability theory. For example, it has been proved |GM05[ Th. 2] that 
if the spectral law (i.e. uniform distributions on eigenvalues) of the self-adjoint matrix 
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An converges to a compactly supported probability measure fj, and Fn = diag(l, 0, . . . , 0), 
then for 6 small enough, 

29 

(1) -logE{exp[n9TiiFnUnAnU:)]} ^ ^ / ' R^m^, 

where i?^ is the so-called R-transform of /i and /3 = 1 or 2 according to wether we 
consider real or complex matrices. The /^-transform is an integral transform of probability 
measures on M. Its main property is that it linearizes the additive free convolution ffl, 
the binary operation on probability measures on M which can be defined by the fact 
that for A, B large self-adjoint random matrices with spectral laws tending to fi^, fiB (as 
the dimension goes to infinity) and U a Haar-distributed orthogonal or unitary matrix 
independent of A and B, the spectral law of A + UBU* tends to /i^ H fJ-B'- the free 
convolution ffl can be thought as the analogue, for the spectral laws of certain large 
random matrices, of the classical convolution for real random variables. For all probability 
measures /x, on M, we have 

(2) RtMvif) = R^l(t) + Rvit) (for if: in a neighborhood of zero). 

Hence in the analogy between the free convolution and the classical one, the i?-transform 
plays the role of the log-Laplace transform, and ([T]) gives a concrete sense to this: the 
/^-transform (more specifically its primitive, which also satisfies ([2])), is the limit of a 
certain sequence of log-Laplace transforms. 

Let us now describe the content of our paper. For each A G [0,1], another free con- 
volution, denoted by Ha and called the rectangular free convolution with ratio A, does 
the same job as ffl for the singular laws (i.e. uniform distributions on singular values) of 
rectangular n x m random matrices which dimensions n, m tend to infinity in such a way 
that n/m tends to A: for n, m large integers such that n/m ~ A, for A, B some nxm real 
or complex matrices with singular laws u^, vb and f/, V Haar-distributed orthogonal or 
unitary matrices independent of A and B, the singular law oi A + UBV is approximately 
va Ha t-'b (see |B09aj or the introduction of |B09b] for a more precise definition of Ba)- 
Like the /^-transform for ffl and the log-Laplace transform for the classical convolution, 
an integral transform linearizes Ha- It is called the rectangular R-transform with ratio A 
and is denoted by (7*^^^: for all probability measures /x, u on [0, -f-oo), we have 

(3) C'iS,^(^) = C\^\t) + Cl^\t) (for t in a neighborhood of zero). 

The main result of the paper gives an interpretation of the rectangular /^-transform (more 
specifically its primitive, which also satisfies ([3])) as the limit of a sequence of log-Laplace 
transforms: we prove that if the singular law of M„ tends to a probability measure and 
n/m tends to a limit A G [0, 1] as ra, m tend to infinity, then for 9 small enough, for a 
sequence of m„ x n elementary matrices, 

(4) -logE{exp[vW^3fJ(Tr(E„?7„M„Kn))]} ^ /? / ' dt. 

11 n,m-s>oo Jq t 

Let us mention that free probability theory has initially been built in the area of operator 
algebras and that concrete relations between free and classical probability theory, like the 
ones of ([1]) and (jl]), are not that common. 



RECTANGULAR i?-TRANSFORM AND RECTANGULAR SPHERICAL INTEGRALS 



3 



Let us also mention that expectations of the exponential of traces of polynomials of 
constant matrices and uniform orthogonal random matrices, which have been extensively 
studied in physics and also other areas, like information theory, are often called spherical 
integrals. See e.g. |Z97t IGZ02t [G09j and the references above for the case of square 
matrices and |SW03t IGTOSj for the case of rectangular matrices. 

The paper is organized as follows. In Section [H we state the main result of the paper. 
Theorem II. 2 [ and discuss it. In Section [21 we recall the precise definition of the rectangular 
i?-transform and prove a result of continuity of the map (A,/i) i — Cl^\ At last. Section 
ini is devoted to the proof of Theorem 11.21 following the ideas of the proof of |GM05t Th. 
2]. 



1. Main result 



1.1. Statement. Let us consider, for all n > 1, an integer m„ > n such that, as n tends 
to infinity, n/rrin tends to a limit A G [0, 1] and an n x real or complex nonrandom 
matrix M„ whose singular values are strictly bounded, uniformly in n, by a constant 
K and such that, as n tends to infinity, the singular law of M„ converges weakly to a 
probability measure that we shall denote by /i. Let us define, for G M, 

(5) In{e) = -logE{exp[y^^3fJ(Tr(E„?7„M„V;))]}, 

n 

where Vn are independent Haar-distributed orthogonal or unitary (according to whether 
M„ is real or complex) matrices with respective sizes n x n, nin x m„ and denotes an 
mn X n elementary matrix (i.e. a matrix which entries are all zero, except one of them, 
which is equal to one). 

In the case where A = 0, we also suppose that there is a < 2 such that 

(6) for n large enough, m„ < 



Remark 1.1. Let K be either M or C according to wether we consider real or complex 
matrices. In{d) can also be considered as the Laplace transform of a certain scalar product 
estimated at a pair of independent random vectors, one of them being a uniform random 
vector of the unit sphere of and the other one being the projection, on K", of a uniform 
random vector of the unit sphere of K™" . Indeed, let us denote the singular values of M„ 
by fin,i, ■ ■ ■ ,/in,n and introduce (see |HJ85j ) some orthogonal or unitary matrices Pn,Qn 
with respective sizes n x n, nin x m„ such that 

>n,i ••• 0" 



Qn- 



fJ'n,n ■ ■ ■ 0_ 

Let also, for each n, {in,jn) be the index of the non-null entry of En- Then the jnth row 
(resp. i„th column) Un = (m„,i, . . . , m„,„) (resp. Vn = (t^n,i, • • • , f„,m„)*) of f/„P„ (resp. 
QnVn) is uniformly distributed on the unit sphere of K" (resp. K™") and one has 

I " 

(7) IniO) = - log E{exp [y/nninO ^("^ Un,kfJ'n,kVn,k)]} ■ 

^ k=l 
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The main result of the article is the following one. 

Theorem 1.2. Set /3 = 1 or /3 = 2 according to wether we consider real or complex ma- 
trices. The function In converges uniformly on every compact subset of {—(3K~^,/3K^^) 
to the function 

Jo 1^ 

where denotes the rectangular R-transform of jj, with ratio A (its definition is recalled 
in SectionlE below). 

Remark 1.3. Note that the function C^^^ is analytic on {—K~'^,K~'^) and vanishes at 
zero, so / is actually well defined and analytic on {—K^^,K~^). 

1.2. Particular cases where the matrices M„ are square (A = 1) or asymptoti- 
cally fiat (A = 0). Let us recall that the R-transform of a probability measure u is the 
function 

Rjz) = G-Hz) - for GJz' ^ "^^^^^ 



z J z — t 

(the convention we use is the one of the analytic approach to freeness |HP00t IAGZ09] . 
which is not exactly the one ofthe combinatorial approach |NS06] : ^combinatories ^^^^ _ 

Let fig be the symmetrization of fi, defined by yUs(^) = for all Borel subset 

A of M, and /i^ be the push- forward of by the function t \^ t^. 

Corollary 1.4. In the particular case where A = 1 (resp. \ = 0), the limit I of In can be 
expressed via the R-transform of fig (resp. fi"^) in the following way 

I{d) = 13 r R^Mdt, (resp. I{d) = 13 r tR^2{t^)dt.) 



Proof. It suffices to prove that C ^^^(t^) = tR^^{t) and that Gj^\t) = tR^2{t). The second 
equation can be found in |B09a[ Lem. 3.2 or Sect. 3.6]. The first equation follows from 
the fact that for all A, Cl^^ = C^^^ and from the fact that for all symmetric probability 
measure z/, by p09al Sect. 3.6], cP{z^) = zR^{z). □ 

1.3. Possible extensions of Theorem 11.21 

1.3.1. Cumulants point of view. For A G [0,1], the rectangular free cumulants with ratio 
A of /i have been defined in |B09at Sect. 3.4] (see also |B07at Sect. 2.2]): this is the 
sequence {c2k{^J'))k>l linked to the moments of /i by |B07b| Eq. (4.1)]. Recall also that for 
X a bounded real random variable, the classical cumulants of X are the numbers Clk{X) 
defined by the formula 

logE(e-) = 5:^z^ 

fc>i 
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Differentiating formally the convergence Inid) — > would get the following 

"classical cumulants interpretation" of the rectangular free cumulants with ratio A: for 
all positive integers k, 

(8) C2k{fi) = lim (^"^")' /^'' I ChkmMEnUnMM))- 

n-s>oo n [2k — 1)1 

This formula can be considered as a "rectangular analogue" of [C03[ Th. 4.7]. According 
to the author's intuition, there could be two ways to prove rigorously. One of them 
would be to use the Weingarten calculus, developed by Collins and Sniady, for the com- 
putation of the expectation of moments of the entries of the matrices f/„, Vn (as in the 
proof of |C03l Th. 4.7]). The other one would rely on the extension of Theorem 11.21 to 
complex values of 6 and notice that the functions in question there are analytic in 6 (so 
that their convergence implies the one of their derivatives). 

1.3.2. Case where M is chosen at random. If M„ is also chosen at random, independently 
of Un and Vn, and the expectation, in ([2]), is taken with respect to the randomness of Un, Vn 
and Mn, then Theorem 11.21 stays true in certain cases (for example if M is a standard 
Gaussian matrix divided by y/nin) , but it can easily be seen that Theorem 11.21 is not true 
in general anymore. However, if the expectation, in ([5]), only concerns the randomness of 
Un and Vn, then /„(6') is a random variable, and for certain sequences of random matrices 
Mn, more than Theorem 11.21 can be said about its convergence. An important example is 
given by the case where M„ = An + PnBnQn with An, Bn deterministic matrices having 
limit singular laws [xa^I^b and Pn,Qn Haar-distributed orthogonal or unitary matrices 
with respective sizes n and m„. In this case, one can prove (with technical additional 
hypotheses), that almost surely, 

lim„^oo ^ log / ev^^«(T'-(i^nC/nM„y„))^f/^^y^ _ 

lim„^oo ^ log / ev^^«(Tr(£„c/„A„y„))^f;^^y^ + ^ log / ev^^«('^(^"^"^"^"))df/„dK- 

To prove it, the ideas are the same as the ones of [GM05[ Sect. 6.1 and 6.2]: first prove 
that the random variable 

-log /"ev^^*™"^^('^(^"'^"*^"^"))d[4dK 
n J 

concentrates around its mean (to do it, use [AGZ09t Cor. 4.4.30] instead of |GM05t Lem. 
24] in order to be allowed to consider the complex case) and then a the ideas of [ GMOSt 
Sect. 6.2]. Since the singular law of M„ converges almost surely to /i^ Ha fJ'B, this gives 
a new proof of ([3]). 

1.3.3. Strong continuity property for the rectangular spherical integrals. One other way 
to extend this work, suggested to us by a referee, would be to use fl25]) to prove a strong 
continuity property for the rectangular spherical integrals, in the spirit of Proposition 2.1 
or Lemma 2.3 of [MOT] . 



M,<z)= ^^dMt)=/ --^dMt)-i {ze[o,K-')). 
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2. Preliminaries about the rectangular _R-transform 

Let yU be a probability measure on the real line which support is contained in [—K, K], 
with K > (we do not suppose /x to be symmetric, how it was the case in the initial 
definition of the rectangular /^-transform) . Let us define the generating function of the 
moments of fi"^ 

t^z , , , f 1 

1 - t^Z JtgR 1 - t^Z 

It can easily be proved that M^2 is nonnegative and non decreasing on [0,K~'^). Let us 
define, for A G [0, 1], TW(^) = {Xz + l)(;z + 1), and 

Hl'\z) = zT^'\M,.{z)). 

Then H^^^ defines an increasing analytic diffeomorphism (in this paper, for I an interval, 
we shall call an analytic function on I a function on I which extends analytically to 
an open subset of C containing /) from [0, onto the (possibly unbounded) interval 

[0,limz^K-^ hI^\z)) such that 

i/W(0)=0, d,Hl'\0) = l, Hl'\z)>z, \im H^^\z) > 

We denote its inverse by h'"^^ . Moreover, T^^^ defines an analytic increasing diffeomor- 
phism from [— l,+oo) to [0,+oo), thus one can define the rectangular R-transform with 
ratio A of /x: 

(9) g^H^)=TW'M ' I for ^7^0, and ^^(0) = 0, 

which is analytic and non negative on the interval [0, Mm ^ Hj^^\z)) (which always 
contains [0, K~^)). 

By Theorems 3.8 and 3.12 of |B09a] . the rectangular /^-transform characterizes sym- 
metric measures, and for all pair /xi,/i2 of compactly supported symmetric probabihty 
measures, fii ffl^ is characterized by the fact that in a neighborhood of zero. 

The following theorem states the continuity of the mapping (A,/i) i — > Cl^'^ in a way 
which is quite different from the one of Theorem 3.11 of [B09aJ (where A was fixed). 

Theorem 2.1. Fix K > 0, let fin be a sequence of probability measures on [—K,K] which 
converges weakly to a limit fi, and let A^ be a sequence of elements of [0, 1] which converges 
to a limit A G [0, 1]. Then the sequence of functions C^^"'' converges to Cj^"^ uniformly on 
every compact subset of [0, K^"^). 

Proof. Recall that C\^^ is defined by OH]). Since, by Heine's Theorem, (A, z) (-)■ T'^'^^ ^{z) 
is uniformly continuous on every compact subset of [0,1] x [0,-|-cx3), it suffices to prove 
that f_i — converges to — uniformly on every compact subset of [0,/^"^). 
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Claim a : For each compact subset E of 'C\[K~'^ ,-\-oo), there is a constant such that 
for any law v on [—K, K], for any c e [0, 1], for any z & E 

\T^'^\M,.{z))\<kE. 
Indeed, for z G C\[i^~^, +oo), for any law v on [—K, K], for any c G [0, 1], 

r(^)(M..(z)) = / TrKMmcv + (1 - c)5,){t'l 

J{t,t')ehK,Ky^ (1 - zt^){^ - zt' ) 

thus kE = max{|l — zt^\~'^ ; \t\ < K^z E E} is convenient. 

Claim b : As v varies in the set of laws on [—K,K] and c varies in [0,1], the set of 
functions 

ze[0,K-')^ ' 

Hi'^ (z) 

is relatively compact for the topology of uniform convergence on every compact subset of 
[0,K~'^). By Ascoli's Theorem, to prove Claim b, it suffices to prove that this family is 
uniformly bounded and uniformly Lipschitz on every compact subset of [0,K~'^). Let us 
fix i/ a law on [—K, K] and c e [0, 1]. Note that we have 



Since, moreover, for all z G [0, Hjf^ (z) < z (indeed, for all z G [0, -ftT"^), Hu\z) > 

z) , it suffices to verify that the sets of functions 

{z ^ ; ly law on [-K, K], c e [0, 1]} 

and {z ^ '"^Jp"!^^^ ; ^ law on [-K, X], c G [0, 1]} 

are uniformly bounded on every compact subset of [0,i^'~^). The family of functions 

ihLM = T(^)(M^2(z)), indexed by z/, c, is a family of analytic functions on C\[-ft' ^, +oo) 
which is uniformly bounded on every compact subset of CWK^"^, +oo) (by Claim a). As 

a consequence, the family of the derivatives " ^ is also uniformly bounded on every 
compact subset of C\[i^~^, +00). Since 

zHjf^'{z)-Hi'\z) _ 1 



and hI'^^ (z) > 1 on [0, X ^), Claim b is proved. 

Hence one can suppose that — converges to a function / uniformly on every 

compact of [0,K~^). Let us fix z G [0,/^"^) and let us prove that f{z) = — (jpr — • If 

2; = 0, it is clear (since all these functions are implicitly defined to map to 1). Suppose 
that 2; > 0. Note that f(z) ^ 0, because for all n, f_i — > 1. Let us denote I — 
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It suffices to prove that / = H^^ (z), i.e. that Hj^'\l) = z. Since 

= hm Hl'\Hl':r\,^^^ 

it suffices to prove that converges to H^^^^ uniformly on every compact subset of 

[0,i^'~^). But it is easy to see, using |AGZ09l Th. C.ll], that M^2 converges to M^2 

uniformly on every compact subset of [0, and then that H^^""^ converges 

uniformly on every compact subset of [0,/^"^). □ 



3. Proof of Theorem 11.21 

3.1. Preliminaries. We shall use the following lemmas several times in the paper. Let 
II ■ II denote the canonical euclidian norm on each Mf'. 

Lemma 3.1. Let (Gj)j>i he a family of independent real random variables with standard 
Gaussian law. Let T he fixed and let, for each n, (o"„,i, . . . , (Tn,n) £ [0, T]" he such that 



(10) ^E<' = i- 



n 

i=l 



Let us define, for each n, X„ = ((Xn^iGi, . . . , cr„^„G„). Then for all k G (0, |), for n large 
enough, 

\Xn\\ -Vn\> n^""} < 2T^"2"-^ 



n 



Proof. Note that by fflOl) . the random variable Nn := 1 is centered. Moreover, 

Var(A^„) = -^-7^^ ^"=1 cr^ i < It follows, by Tchebichev's inequality, that for all 

«:e(o,i), 

P{|iVn| > n-""} < 2TV'^-\ 
To deduce that for n large enough, 

\X„\\ 



P 



n 



it suffices to notice that the function yf- is 1-Lipschitz on [1/4, +00) and that n < 3/4 
for n large enough. □ 

Lemma 3.2. Let ji he a prohahility measure whose support is contained in [—K,K], fix 
A G [0, 1], ^ G [0, K-^) and define 7 = d^^\e^). Then 

(11) M,. f = 7- 



TW(7) 

Proof. By the definition of Cj^^^ given in (xy-^^ 2) ~ ^^'^^(t)' hence ^^f)^^^^ = 

H^^^ (6*^). Since 7 > 0, rp(\^^-^ £ [0, -ft'^^) and one can apply the function hI^^ on 
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both sides. We get i/^^^ (r(%y) = i-^- 



It follows that T(^) (^M^2 (^t(S^)) = ^^^H^)- Since both M^2 (5xxr(^) and 7 are non- 
negative real numbers, one gets fill I) . □ 

Lemma 3.3. Let Xn be a sequence of nonnegative random variables, with positive expec- 
tations. Let Zn be a sequence of real random variables such that there exists deterministic 
constants C, 77 > such that for all n, \Zn\ < Cn^~^. Then as n tends to infinity, 

-logE(X„e^") = -logE(X„) + 0(1). 
n n 

Proof. It suffices to notice that we have Xne"*^"^"" < X„e^" < X^e"^"-^'". □ 



3.2. Notation for the proof of Theorem 11.21 In the next sections, o(l) shall denote 
any sequence of functions on (— i^^^, K^^) which converges to zero as n tends to infinity, 
uniformly on every compact subset of (— i^"^, K~^). Also, we shall work with the notation 
introduced in Remark II. II and handle /„(6') via Formula ([7]). 



3.3. Case where Af„ has null singular values. Let us suppose the result to be proved 
in the particular case where for all n, /c, iXn,k > 0, and let us prove it in the general case. 
We set, for each n, k, 

P'n,k = fJ^n^k + min{m^^, {K - /i„,fc)/2} 
and define the perturbation of In{d) 

ln{0) = - logE{exp[v/nm„6' ^{Un,kP'n,kVn,k)]} ■ 

k=l 

The uniform law on the jln,k^ converges weakly to /i as n tends to infinity and we have 
< jln,k < K, so by hypothesis, if follows that lim„^oo ln{0) = H^)- Note moreover that 

I " 

Inid) = -logE{exp[v/ri^6'^3fJ(M„,fe/i„,fct'„,fc)]e^"}, 

k=l 

with 

n 

Zn = ^/nninO ^ 3ft(n„,fc^^n,fc) min{m^^ {K - /i„_fc)/2}. 

k=l 

Since |Z„| < K~^, Lemma [3.31 allows us to claim that lim„_!.oo /n(^) = 
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3.4. Deducing the complex case from the real one. Let us explain how one can 
deduce the complex case from the real one. Let us use, in this paragraph, the notation 
In,rnn{(^, Mn) to cmphasizc on the value of each of the parameters defining /„(^). We have 



- log E{exp[y/nmn9 ^ ^iUn,kfJ'n,kVn,k)] } 

^ k=l 
I " 

- log E{exp[v'nm„6' /i„,fc(3?(Mn,fe)3^(^^n,fc) - S(^in,fc)5>(^n,fc))]} 



(12) 



n 



21 



k=l 



(1) 

2n,2m„ 


















Indeed, for {zn^i, . . . , Zn,n) a vector with uniform distribution on the unit sphere of the 
vector (3fJ(zi), . . . , 3fJ(z„), 9(zri)) has uniform distribution on the unit sphere of M^" 

(use the realization of such vectors as the projection of standard Gaussian vectors on the 
sphere to see it). 

"M„ 
-M„ 



Since the singular values of 



are the one of M„ with multiplicity multiplied 



by two, (fT2|) allows to deduce the complex case from the real one. 



Let us now prove Theorem 11.21 in the case where all ^n,k^ are positive and where /5 = 1. 



3.5. Proof of Theorem II. 2t a) Expression of /„(6') with a Gaussian integral. As 

written above, we suppose from now on that all /i„,fc's are positive and that /3 = 1. In 
this section, we shall first explain how to replace, in Formula ([7]), \/nUn,k and ^Jm^Vn^k by 
independent standard Gaussian variables (Formula f lT7|) ) and then inject (quite artificially 
first) Ci^V^) in the formula of 4(^) (Equation ([21])). 
For each n, let us define the function 

n 

fn : ((a;i,...,x„),(?/i,...,?/^J) G M" x M™" i — > ^[in^kXkVk- 

k=l 

By |GM05l Sect. 1.2.1], up to a change of the probability space which does not change 
the expectation, one can suppose that there are independent standard Gaussian vectors 
Xn,yn of M", respectively M."^", such that 

-^n Vn 

Un Tj TT) Tj 77- 

1 1 "^n II II Vn 1 1 

Let us fix K G (0, 1/2). If A > 0, the precise choice of k G (0, 1/2) is irrelevant, but if 
A = 0, we choose k G (^^, |) (a is the one of (^). Let us now define the set 

An := \{x,y) G R" X M™" ; |||x|| - y/n\ < n^~'^ , \ \\y\\ - ^Jm^\ < m-F'^j . 
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The event {{xn,yn) £ ^n} is independent of {un,Vn) (because the density of the law of a 
standard Gaussian vector is a radial function), thus 

ln{0) = ;^logE [l{x„,y„)eAnexp{y/nmJfn{un,Vn))] - ^log 
Moreover, by Lemma [STTl P{(x„,y„) G An} — > 1 as n — oo, thus 

(13) ln{0) = -logE [l{x„,y„)eAn exp{y/nmn9 fn{Un, Vn))] + o{l) . 

Moreover, note that on the event {(a;„,?/„) G An}, 



^rn^-m^ " < \\yn\\ < y/m^ + 



thus, since > n, 



(14) ^JmK^ — ?)^frrQn'i < ||y„,|| < ^nmn + 3A/m„n2 

If A > 0, since mn/n is bounded, it follows that there is a deterministic constant C 
independent of n such that on on the event {(a;„, ?/„) G v4„}, 

(15) I IknII llZ/nll - V^"^ri | < Cn^~'^ . 

If A = 0, it follows from fll4p and ^ that for rj = + /t (which is positive by definition 
of k), for n large enough, 

(16) I \\yn\\ - y/nnin I < 3n^"''. 
Note that by 

Ud) = ilogE[l(.„,,„),^„e^^"(^-^")+'fe^(^-"^'"""^"")] + 0(1), 

< K. Hence by Lemma 



II II II ?/n II 



and that for all n, A;, < ^k,n < which implies that 
ESland (USD (or ([16]) if A = 0), 

(17) UO) = ^logE [l(.„,,„)eA„e^^"(-'^")] + o(l). 
Note that on the event {(x„,?/„) G ^n}? we have 

- 272^"" < n - 2nm~'' + nm"^'' < ;^||y„|P < n + 2nm';^'' + nm:^'^'^ < n + ?,n^~ 
thus for all n, on the event {(a;„,?/„) G 

(18) ||x„f - n + ||?/„f - n < 6n^"'". 

Now, let us define, for each n, 

(19) ln{e) = cl:-\9') for /i„ = -5^<5,„.,. 

fc=i 
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Note that /i„ is the singular law of M„, which tends to /i. Hence by Theorem 12. ![ we have 
(20) 7„(e) ^ Cl^\0^) uniformly on every compact subset of {~K 

SO by (IT5]1 . for every such compact set E, there is a constant Q^; such that for all n, for 
all 9 E E, on the event {(a;„,|/„) G ^n}, we have 



n 



hni6)i-\\XnV + 



n) \ < QeU 



1-K 



Hence, by ( 1T7|1 and Lemma |33 



/„(0)=-logE 



l(xn,j/n)eA„exp <j 0/„(x„,?/„) -7n(^) ( 2ll^"ll + 2r7r"^"' 



=7„(^) + -logE 
n 



+0(1). 



V 

■■ = Jr^{e) 



Thus, by (|20 
(21) 



in{e) = c^^\e') + -\ogUd) + o{i) 

^ n 



3.6. Proof of Theorem II. 2t b) Convergence of the Gaussian integral. We have, 
assimilating the vectors of M" and M™'" with column-matrices. 



(22) Ud) = (27r)-^ 



for 



lA„(a;,i/)exp<^ -- [x* y*] T„ 



dxdy. 



an{0)ln An{9) 
T„ := An{9) bn{9)In On,m„-n 

Om,„— ra,n Om„— n,ra bn(^9)Im„—n 

where a„(6') = 1 + 7n(^), &n(^) = 1 + ^7n(6') and A„(6') is the diagonal nx n matrix with 
diagonal entries 

Notation: In this section, in order to lighten the notation, we shall write J„ for Jn{9), 
an for an{9), etc. We shall also use the notation of the matricial functional calculus (thus 
assimilate a and a/„, etc.). 

Lemma 3.4. Let us fix n > 1 and let a, b be real numbers and A an invertible diagonal 
real n x n matrix. Let us define 



1 



2A±2{b-a)y/A 



and 



T 



a A 
A b 



D 



P 



2A/+ 2A/- 
_(6-a)/+ + v^/+ (6-a)r-VA/- 

Then P is an orthogonal matrix and we have T = PDP*. 



r+ 
r- 
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Q 



H 



/+ 

/- 



Proof. One can easily verify that P is orthogonal. Let us define 

2A 2A 

{b-a) + ^/A {b-a)-VA_ 

Then P = QH. One can easily verify that TQ = QD. If follows that TQH = QDH. 
Since HD = DH, TQH = QHD, i.e. TP = PD, thus T = PDP\ □ 

For 6* 7^ 0, let us define A„, r^, as in the lemma, using A„ instead of A, a„ instead of 
a and 6„ instead of h. Let us define P„ in the same way, extended to an (n + m„) x (n+m„) 
matrix by adding Im„-n on the lower-right corner, i.e. 



2A„/+ 2AJ- 
















and Dn extended to an (n + m„) x (n + m„) matrix by adding bnlm„-n on the lower-right 



corner, i.e. 



r+ 
r- 
hj„ 







For 6^ = 0, we set = 1, P„ = = 4+m„- 

Let us denote, for X an (n + m„) x (n + m„) matrix, X{An) = {X 



{x,y) e An}. 



Let us also introduce a standard Gaussian random column vector in M"+'^" ^ that we shall 
denote by 



'7^ 7+ 7- 7- 7O 7 

^■^n,l5 • • • ) ^n,m -^n,!' • • • 5 ^n,ni ^riA^ . . . ^ ^, 



y 

n,m„—n) " 



— '^n 



=Z0 



We have, by ([22]) and Lemma [331 

J„ = (27r)-^ f exp{-^ [x* y*] P^D^P, 

J An 

Thus, since Pn is an orthogonal matrix, 



n + mn 
2 



[ exp{-^ [x* y*] P>„ 



}dxdy. 



}da;d|/. 



Jn = (27r) 
Hence, by definition of Dn, we have 

n 

Jn = (2vr)-^[e--n^ 
which, by definition of Zn, can be written 

n 

(23) Jn = [C""" n(^"^" - P{Z„, G /d;P^(A„)} 



71, r n,ij 



/ exp{-^(||s||2 + llylHldxd?/, 
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Let Xn = . . . be the vector of the first n coordinates of PnDn^^'^ and 

Yn = (^,1, • ■ ■ 5 Yn,m.„y be the one of the m„ last ones. By definition of the set An, we 
have 

(24) PnD-^/^Zn eAn ^ \\\Xn\\ " < n"'' and |||K|| - < m-\ 

Claim : Both events {\\\Xn\\ — \/n\ < n^'^} and {|||l^n|| — < m^^j have prob- 

abilities tending to one as n tends to infinity, uniformly on every compact subset of 
{—K~^,K^^) (the random vectors X„ and Yn depend indeed on 6). 

For = 0, Xn (resp. Yn) is a standard Gaussian random vector of M" (resp. M™'" ). For 
7^ 0, by the definitions of P„ and Dn, 

2An(/„+(r+)-^/^Z+ + f-{r-)-'/'Z-), 

-{bn -an + v^)/+(r+)-V2z+ + {bn -an' VA^)/- (r")" V^Z" 

b-'^'z2 



Xn 

Y 

J- T). 



Thus for each n, Xn (resp. Yn) has the law of 

(0"„,lGl, . . . , CTn.nCn) (resp. ((T^ ^^Gi , . . . , (T^ ,„^Gm„ ) ) , 

for {Gi)i>i a family of independent real random variables with standard Gaussian law and 
where if 6* = 0, all an/s and a'n /s are equal to 1, whereas if 6* 7^ 0, for each i = 1, . . . ,n, 



4AL- 



(2A„,i + 2(6„ - a„)^A„,i)(a„ + 6„ + 

2 



+4A 



(2A„,i - 2{bn - a„)A/A~)(a„ + - a/A~) 



for each i = n + 1, . . . , m„, a'n i' = bn^ and for each i = 1, . . . , n, 



/ 2 



2(6„ - a„ + a/A^ 



2(v^- (&n-an))^ 



(a„ + bn + VKn)[2{bn - an)V^ + '^^n] (^n + &n " VA;;)[-2(&n - a„)VA;:+ 2A„] 

Hence by Lemma [3.11 to prove the claim, it suffices to prove: 
(25) Ve > 0, 



sup 

\e\<K-'^-e l<i<n 
n>l 



(26) 



\/e > 0, 



sup sup a. 

\0\<K-'^-s l<i<m- 
n>l 



sup an i < +00 and - . = 1, 

- n — ^ 

i=l 

n,i < +00 and — V a; / = 



i=l 



Note first that (ES]) and (ES]) both hold when 6 = 0. 

A straightforward computation leads, for 7^ 0, to the formula 
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But (removing the indices) 

{b^ + A-ay -Ab^A = {2b^ -2ab + AX^y -Ab\b^ -2ab + a^ + AX^) 

= 16A^-16a6A2. 

It follows, writing 7„ for 7n(6'), that 



an&n-A^^, 1 + 7. T(^H7n)-^^X. 1 + ^" ^-rf^'"'''' 

By definition of jniO), we have 7n(6') > 0, hence T^"^\'jn{d)) > 1- Since for all n,i, 
|/^n,i| < it follows that the first part of ( 125|) holds. Moreover, by the definition of fin 
given in (fTOl) . we have 



1 " 1 

l + 7r 



9^ 



-M ' 

By (fTTl) . it follows that the second part of (l25l) also holds. 

Let us now prove (126|) . A straightforward computation leads, for 6* 7^ and i < n, to 
the formula 

/ 2 O-n 



O'nbn — A, 



2 • 



Hence for 7^ and i < n, 



2 1 1 



(tL f = 7^ X 



whereas for i = n + 1, . . . , m„, cr^^ / = ^ — r^y. The first part of fl26|) holds for the same 
reasons as the first part of fl25l) above. Moreover, writing 7„ for 7n(6'), we have 



1 / 2 n 6^ \ n , rrin-n 

o"r, ,- = 7^- — r^u2 „ , , — - + 



i=l 



.(l + ^7n) ""{T^^Hln)) ^"(l + ^7n) m„,(l + ^7n) ' 

By (ITTi) . it follows that the second part of also holds. 

The proof of the claim is complete. As a consequence, by (^^, the probability of the 
event {PnDn Zn E An} tends to one as n tends to infinity, uniformly on every compact 
subset of {—K~^, K~^) (remember indeed that the matrices P„ and Dn depend on 6). So 
by ([23]), we have, still writing 7^ for 7„(6'), 

(28) ^log(J„(0)) = 



n-run , n log(T^'"n^(7„)) If r 21 1 , /1^ 

,^ n log(l + 7„) If 21 1 u\ , n^ 

m„ 2 ^ JtG\-K.K] T^™"M7r,, 
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By hypothesis, which is defined in ( 1T9|) . converges weakly to /i. Using (120|) and 
[AG Z 09 1 Th. C.ll], one easily sees that, writing 7 for Cl^\6^), we have 

- log( J„(^)) = log(l + A7) - Mill) _ 1 /■ log{l - -^t^}d/.(t) + 0(1), 

where in the case where A = 0, 2^1og(l + A7) has to be understood as |. By ( |2T|) . one 
gets, still writing 7 for Cj^'\6^), 

/..(«) = 7 - ^ log(l + A,) - !25(L±2) _ 1 £ leg{l _ ^*nd.(*) +o(l). 



/(O) = /(O) = (indeed, by Q, Q^(0) = 0). So to conclude the proof of Theorem Ol it 
suffices to verify that / and / have the same derivatives on {—K~^, K~^) (equal to 7/6*), 
which can easily be done using (ITT]) again. 

Aknowledgements: The author would like to thank Mylene Mai'da for her comments 
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